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Abstract 

We present a microscopic approach to quantum dissipation and sketch the derivation of the kinetic equation 
describing the evolution of a simple quantum system in interaction with a complex quantum system. A typical 
quantum complex system is modeled by means of parametric banded random matrices coupled to the subsystem 
of interest. We do not assume the weak coupling limit and allow for an independent dynamics of the "reservoir". 
We discuss the reasons for having a new theoretical approach and the new elements introduced by us. The present 
approach incorporates known limits and previous results, but at the same time includes new cases, previously 
never derived on a microscopic level. We briefly discuss the kinetic equation and its solution for a particle in the 
absence of an external field. 
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1. Introduction 

The notion of a complex (quantum) system 
might not be an uniquely and universally ac- 
cepted term, but it is certainly a very widely 
(over)used one. We shall call systems with a large 
number of degrees of freedom complex systems. 
One might reasonable ask: "Why a new term?" . 
The real reason we believe, besides the tempta- 
tion to add something "new" to the jargon, is 
to make a distinction with much simpler systems 
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studied in the near past by the community con- 
tributing to this volume. Typically, in complex 
quantum systems a relatively small number of de- 
grees of freedom (DoF) can be singled out in one 
way or another and one might be interested, or 
forced, to study the dynamics of this reduced set 
of DoF. Since these selected DoF are coupled, ei- 
ther weakly, but most likely strongly, to the rest 
of the DoF of the system and even to the rest of 
the Universe, one is faced with the problem of de- 
vising ways to predict and describe the dynamics 
of non-isolated quantum systems. 

In a certain sense such systems can be also 
called open quantum systems. However, the term 
open quantum system has a slightly different con- 
notation. In the case of an open system there is 
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an "orifice" through which once in a while the 
system can "leak out" , or better put, "drop out" , 
and never come back. In an open system the 
role of the "rest of the Universe" is typically re- 
duced to a simple "black hole" or a "drain" or a 
"sink" , which "swallows" matter out of the sys- 
tem as is "been fed" , thus "stuff simply falls out 
of the system", and the "rest of the Universe" 
does not really exert any real force on the sys- 
tem of interest. One can say that the "rest of 
the Universe" is "forced fed", but with no "real 
effort" . The "rest of the Universe" can also serve 
as a "source" of matter, but again, this is done in 
the same somewhat "non-obtrusive way" as the 
time-reverse process of "swallowing" . The "rest 
of the Universe" can even play both roles, "sink" 
and "source" , at the same time. Thus an open 
system is nothing else, but a system with some- 
how modified boundary conditions. One might 
argue that a matter flux can influence dynam- 
ics. That would be correct if one would consider 
a system of interacting particles, whose number 
can vary in time. The folks typically involved 
with open quantum systems study physics at the 
independent particle level however. 

A complex system in the sense we shall be us- 
ing is however somewhat different. We shall not 
allow (yet) matter exchange, simply for the sake 
of simplicity, but allow only for energy exchange 
between the part of the system we have decided 
(or have been forced to) to focus our attention 
on and the "rest of the Universe" . Both parts of 
the system have their own dynamics, but at the 
same time, the interaction between them could be 
strong enough as to modify the dynamic evolu- 
tion in a qualitative or quantitative manner. The 
qualifier "complex system" refers as a rule either 
to the "rest of the Universe" only or to the entire 
system. The dynamic evolution of the "rest of 
the Universe" is routinely assumed to fall in one 
or another universality class, described in Quan- 
tum Mechanics by some random Hamiltonian ^ 
or by some other Hamiltonian with "complex dy- 
namics" . If the interaction between the two sub- 
systems would have been weak, one would have 
dealt with a situation typically studied in Sta- 
tistical Physics textbooks. The energy exchange 
between the two subsystems could be character- 



ized as in Thermodynamics as either mechanical 
work (thus reversible) or heat exchange (conse- 
quently irreversible). If irreversible processes oc- 
cur, one obviously has entropy production as well. 
Moreover, if one would have concentrated one's 
attention only on the dynamics of the smaller sys- 
tem, the bigger one would most likely have been 
called a thermal bath, and one would have aimed 
at obtaining a kinetic description of the dynam- 
ics of the smaller /simpler system. If the "rest of 
the Universe" is huge, one might reasonable argue 
that the influence of the subsystem on it should 
be negligible and thus could be ignored. Besides 
the fact that that is not always the case of inter- 
est (when the "rest of the Universe" is not huge, 
like in the case of the interaction between the col- 
lective DoF in a nucleus and the rest), there is no 
need to influence all the "rest of the Universe" , 
but only the part lying nearby and if the heat 
conductivity of the "rest of the Universe" is low, 
energy will not be quickly transferred to the out- 
skirts of the "rest of the Universe" . 

There was a significant interest over the years 
in quantum problems similar to the type we have 
just described. Caldeira and Leggett [||, using a 
technology outlined by Feynman and Vernon ^ 
have studied the motion a "heavy" particle in- 
teracting with a reservoir of harmonic oscillators. 
The reservoir of harmonic oscillators is assumed 
to mock a quantum "heat" reservoir. The quan- 
tum Hamiltonian describing the "Universe" ac- 
cording to Caldeira and Leggett has the form 

H = ^ + UiQ) 

tj2.n.+ 2 J' 

with the corresponding p's and g's and so forth 
and the coupling constants Cfc. A continuous dis- 
tribution of these harmonic oscillators leads to 
a dissipative dynamics of the "heavy" system de- 
scribed by momentum P and coordinate Q and in 
particular to a Langevin description in the clas- 
sical limit. The fact that this particular classi- 
cal limit should emerge in the high-temperature 
regime, was one of the requirements imposed on 
the structure of the model Hamiltonian in Eq. 
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(y). This model Hamiltonian suffers however 
from one rather simple inadequacy, which is at 
least intellectually unsatisfactory. Let us assume 
that we would like to describe the dissipative 
motion of a particle in a homogeneous medium. 
There is no external potential U{Q) then and it 
is completely unclear in this case how one can 
define a coupling to the reservoir, as there is ob- 
viously no special point with respect to which one 
can measure Q, since the medium is supposed to 
be homogeneous. One might consider instead a 
coupling through a derivative term, i.e. to the 
momentum of the "heavy" particle. Then one 
would have had to alter the time-reversal prop- 
erties of the Hamiltonian for no other particular 
reason. A coupling to the momentum typically 
implies some gauge field. In particular, such a 
mechanism will not work in one dimension, as 
no magnetic field exists in ID. Moreover, simply 
having to consider a gauge field for the sake of 
having dissipation looks like an overkill, at least 
from a theoretical point of view. Coupling to spin 
degrees of freedom for the same reason should be 
ruled out as generic. Thus we seem to find our- 
selves in a quandary. Even in the presence of an 
external field, the fictitious harmonic oscillators 
become more and more elongated as the "heavy" 
particle moves away from the "center" and thus 
the "reservoir" becomes "infinitely stretched" and 
consequently, due to the presence of an infinite 
number of DoF in the reservoir, "infinitely ex- 
cited" too. Moreover, even though it is entirely 
reasonable to be guided by the "known" classi- 
cal limit when "guessing" the structure of the 
model Hamiltonian, one can hardly accept this 
as a satisfactory recipe and a suitable solution 
for all cases. On one hand, the form of the fric- 
tion force is not unique even in classical physics, 
classical friction force can be proportional to the 
velocity at various powers or even independent of 
velocity. On the other hand, atomic nuclei for 
example represent systems in which there is no 
classical limit and they are certainly not unique 
in this respect. Moreover, there very well could 
be a genuine quantum dissipative regime, which 
has no counterpart in the classical limit. 

We would like to start from a model "typical 
Hamiltonian" of a "typical quantum complex sys- 



tem" , with a minimum and relevant amount of 
input and find out under what conditions a dis- 
sipative dynamics emerges, what is its character 
and how various properties of this dynamics de- 
pend on the particular properties of the initial 
Hamiltonian. It is not obvious that one will ar- 
rive at the same type of description as Caldeira 
and Leggett have obtained. As a matter of fact 
we find that choosing this alternative approach 
one sometimes arrives at a dissipative dynamics 
similar to that found by Caldeira and Leggett. 
The surprise however is that dissipative dynam- 
ics with entirely new characteristics emerges as 
well. The particular example of quantum Levy 
flights and fractional kinetics is one remarkable 
case H]. 

There is another direction of inquiry in the the- 
ory of open quantum systems or quantum dissipa- 
tive dynamics, which we shall mention only very 
superficially Q . The main emphasis in these ap- 
proaches is phenomenological, at the level of the 
von Neuman equation and its possible extensions, 
such as to accommodate a Markovian evolution, 
to preserve the positivity of the density matrix 
and to conserve the probability. As one might 
have expected from a pure phenomenological ap- 
proach, there has been no effort to understand 
the nature and the origin of various terms added 
to the von Neuman equation from any underlin- 
ing microscopic dynamics, as no such derivation 
has ever been attempted, except cases when per- 
turbation theory could be implemented and one 
ends up with a textbook master equation. 

This sketchy review of the literature is not 
meant to be exhaustive or comprehensive in any 
sense. The aim was merely to point the reader 
to some literature leads and to justify in some 
very limited way our approach and line of reason- 
ing. We deeply apologize for not devoting more 
space and emphasis to so much fine work, not 
mentioned here mostly due to space limitations. 
The rest of the contributors to this volume will 
surely complement us. 

2. A generic quantum complex system 

We shall describe here what we view as a 
generic quantum complex system and the rea- 
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sons why we have chosen this perspective. The 
Hamihonian of the "entire Universe" we choose 
as follows Pi 



n=^^+U{Q) + H{Q), 



(2) 



where H{Q) is supposed to describe the "rest 
of the Universe" and the interaction with the 
selected DoF with canonical coordinates {P,Q), 
which could stand for one or more pairs of canon- 
ical variables. The "rest of the Universe" is as- 
sumed typically to be a very large system, with a 
large heat capacity. It is convenient to choose an 
average level density of states of the "rest of the 
Universe" which is locally exponentially increas- 
ing, 



p{s) = Tr<5(H(Q) - e) = po exp(/3e). 



(3) 



The overline stands for averaging over some por- 
tion of the spectrum or for an ensemble averag- 
ing to be described shortly. Thus P = 1/T is 
nothing else but the inverse thermodynamic tem- 
perature of the "rest of the Universe". So far 
we have always chosen an average level density 
which is independent of the coordinate Q 
1^. In physical terms this means that only heat 
exchange is allowed so far between the selected 
subsystem and the "rest of the Universe" . We do 
not see any difficulties of any kind to incorporate 
the possibility for the system to perform work on 
the "rest of the Universe" as well (or vice versa) . 
Allowing for the temperature T to vary appar- 
ently could require some further small technical 
developments of our formalism. The relevant "de- 
pendence" of the Hamiltonian H{Q) on the sub- 
system coordinate Q is implemented by modeling 
it with a random matrix and specifying the 
second cumulants of its matrix elements; H 



[HiX)U[HiY)]^,^-[H{X)],, [H{Y)i 



27rv/p(efe)/o(£i) 



X exp 



2kq 



G 



X-Y 

Xn 



(4) 



Here is the so called spreading width, kq de- 
fines the handedness of the random matrix and 



it can be interpreted as the energy of a "single 
kick" . To Xq we refer to as the correlation length. 
Even though we have used GOE-type of induced 
spectral fluctuations, similar results are obtained 
for the GUE or GSE case. This in a way could 
have been expected. It is natural to expect that 
most/many complex systems have universal spec- 
tral fluctuations at some small energy scales as 
both theory and experiment indicate. A GOE 
Hamiltonian is suitable for an even number of 
fermions, while a GSE Hamiltonian for an odd 
number of fermions . A "reservoir" with an odd 
or even number of fermions could not influence in 
qualitatively different ways a given subsystem. 

The correlation length Xq is perhaps one of 
the most interesting parameters of the entire the- 
ory. The function G{{X - Y)/Xo) (G(0) = 1 
and G{x) = G{—x) = G*{x) < 1) is assumed to 
have a bell shape with a characteristic width Xq. 
Thus Xq specifies by how much one has to move 
the subsystem (Q, P) in order to induce some no- 
ticeable changes in the "rest of the Universe". 
Obviously, if the "rest of the Universe" does not 
react in any way and does not change its proper- 
ties while the subsystem (Q, P) evolves, there is 
no energy exchange and no change in entropy as 
well. (We trust that this is neither the place nor 
the time where one should discuss at any length 
why entropy can vary, even though one might 
claim that the entropy of the entire system is 
strictly conserved.) The restriction adopted here 
of "translation invariance" has no fundamental 
consequences and apparently can be lifted with- 
out much effort. 

From the characterization given so far of the 
Hamiltonian H{Q) of the "rest of the Universe" 
and of the coupling to the subsystem (Q, P) one 
can see the differences with other approaches 
mentioned in the introduction, in particular with 
the most "microscopic" theory of Caldeira and 
Leggett The "translation invariance" of the 
present approach includes the case of the motion 
of a subsystem {Q, P) in a homogeneous medium. 

From this point on the development of the the- 
ory follows Fcynman and Vernon strategy , one 
writes down a double path-integral for the den- 
sity matrix of the entire system, one integrates 
out the "rest of the Universe" and in the end one 
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arrives at an evolution equation for the density 
matrix of the subsystem {Q,P). It is surpris- 
ing that in the hmit of high temperatures and 
large bandwidth of the matrix [H(Q)]ki (namely, 
when Kopis) 2> 1) one arrives and an extremely 
simple looking evolution equation for the reduced 
density matrix of the subsystem p{X^ Y, t), which 
reads: 



dt \2M 2M ^ ^ 



G' 



G 



X 



X -Y 
Y 



X 



Py 







1 



p{X,Y,t). 



U{Y) 



(5) 



Hopefully no confusion with the notation used 
earlier for the average level density should arise. 
The derivation, which is somewhat lengthy and 
requires an extensive discussion of various steps, 
is described in our earlier works and shall not be 
repeated here. The "rest of the Universe" is al- 
lowed to have its own dynamic evolution and it 
is not assumed to be in thermal equilibrium or to 
have a prescribed evolution. Some rather inno- 
cent assumptions are also made about the nature 
of the initial conditions . In the case when 

one simply drives the "rest of the Universe" one 
can show that the energy distribution has signif- 
icant non-Maxwellian distortions The cou- 
pling between the subsystem (Q, P) and the "rest 
of the Universe" is not assumed to be weak at any 
stage, even though the limit of weak coupling is 
obviously embodied as well in the evolution equa- 
tion (|). 

By means of a Wigner transform and in the 
case when for small arguments 



G{x) K. 1 - + . . 



(6) 



the evolution equation (|^) reduces in the limit 
?i — >• to the Kramers equation ^ , and thus the 
equivalent Langevin description is also obtained. 
The classical friction force in this case is propor- 
tional to the first power of the velocity, the fric- 
tion and diffusion coefficients are given by 
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2MXi ' 



(7) 



QQ - -^2^ (8) 

and the standard Einstein dissipation-fluctuation 
relation is satisfied. 

One might have considered the fact that this 
classical limit is exactly reproduced under ap- 
parently correct conditions as a necessary self- 
consistency check of the entire approach. The 
power of the present approach is demonstrated 
however by being able to generate other distinct 
and consistent classical limits as well, which, as 
far as we can judge, was not shown to be possible 
within the framework of approaches previously 
known in the literature. In particular, a rather 
innocently looking change of the correlator G{x) 
for small arguments to 



G{x) = 1 - ixr + 



(9) 



with < a < 2 leads to apparently the first 
microscopic derivation of fractional kinetics with 
Levy flights |Q] , without assuming or introducing 
any kind of exotic noise. The only "noise" is of 
the usual rather boring Gaussian type. A com- 
pletely different behavior, turbulent-like diffusion 





(10) 



is also in the realm of this type of evolution equa- 
tion (||). The double angle brackets stand for cu- 
mulants. 

Surprisingly, in many cases the evolution equa- 
tion can be solved in quadratures for arbitrary 
type of coupling between the two subsystems [^,|[ 
0. Time-dependent solutions of the Schrodinger 
equation are known for the case of free mo- 
tion, linear potential, harmonic potential. Path- 
integrals also can be calculated in cases when the 
action is essentially quadratic in p's and q's. To 
our surprise the solution of Eq.(||) can be ob- 
tained easily if the classical trajectories in the 
absence of dissipation (i.e. no G{{X — Y) / Xq)) 
are known. Several cases were discussed by us in 
Refs. [^,01 . Here we shall illustrate a few salient 
points of our approach on the example of motion 
in the absence of any external field (U{Q) = 0), 
but for arbitrary type of coupling to the "rest of 
the Universe" . In ID and representing the den- 
sity matrix in the coordinates r = {X + Y)/2, 
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s — X — Y one can show that for an arbitrary 
initial density matrix po F) , for t > one has 



p{r,s,t) = 
dr'dk 



PQ\r ,s 



kt 



ik{r — r') 



hk 



ds'[G{s'/Xo) - 1] 



s-kt/M 



(11) 



This form of the density matrix is especially 
suited to evaluate the coordinate and momentum 
cumulants as functions of time. 

In the long time limit the momentum distri- 
butions reaches an equilibrium, characterized by 
vanishing odd cumulants, while the even cumu- 
lants are given by (n > 2) 

for the particular case of a Gaussian correlator 

G(a;) = exp(-2;V2). (13) 

The second cuniulant has the expected value 

= 2Mr, (14) 

if only for small arguments the correlator behaves 
as in Eq. (^). There is nothing special about the 
Gaussian form of the correlator G{x) used here, 
which was simply chosen for illustrative purposes. 
The parameter which governs the deviation of the 
equilibrium distribution from the one expected 
in the weak coupling limit Maxwellian distribu- 
tion is h/Xo and not as one would have naively 
guessed the "strength" of the couphng TK H/Xq 
could be interpreted as the characteristic momen- 
tum exchanged between the system and the "rest 
of the Universe" in "one completed interaction 
act" . The equilibrium momentum distribution is 
narrower than a pure Maxwellian distribution at 
the same temperature. 

Some information about the coordinate distri- 
bution has been inferred as well The second 
cumulant reveals what one would have expected, 
namely that for t oo 



((Q2)) « 2DQQt, 



(15) 



again, if only for small arguments the correlator 
behaves as in Eq. (^. A qualitative analysis of 



the coordinate distribution seem to indicate the 
presence of "longer tails" than one would have ob- 
tained in the case of a simple Brownian particle. 
That does not seem so surprising now when we 
have learned that under very mild circumstances 
one easily obtains Levy flights behavior and frac- 
tional kinetics j|] or turbulent-like diffusion Q. 
However, a comprehensive analysis of the entire 
coordinate distribution is not achieved yet. 

Analytical results can also be obtained for a 
particle in a linear potential or in a quadratic po- 
tential (both harmonic oscillator and parabolic 
barrier) and qualitatively similar deviations from 
the textbook cases of the weak coupling regime 
are found. The very interesting case of tunnel- 
ing in a double well potential has been studied 
only numerically so far |^J^, along with other 
cases which have some resemblance to under the 
Coulomb barrier collision and tunneling of two 
atomic nuclei. Evidently, our approach is not lim- 
ited to ID as one might have concluded, simply 
because we did not explicitly display spatial sub- 
scripts. 

3. Concluding remarks 

On a closer analysis and in retrospective the 
fact that we have used a random matrix ensem- 
ble, with various adds-on, does not seem to have 
played a fundamental role. If that would have 
been the case we would have seen some distinc- 
tion between GOE, GUE and GSE cases, which 
we did not confirm |1C], even though at one time 
that was expected |12], but later on apparently 
not confirmed as well |l^. The fact that a ran- 
dom matrix approach is not critical to obtaining 
the final answer was confirmed as expected in Ref. 
fl^ in the weak coupling limit. As we have men- 
tioned in the Introduction, it would have been in- 
deed very strange if a GOE-reservoir would lead 
to a different kinetic equation than a GUE- or 
GSE-reservoir. 

However, the particular approach based on the 
use of parametric banded random matrices al- 
lowed us significantly more freedom in choosing 
various forms of coupling between the "rest of 
the Universe" and the subsystem of interest. In 
the wide band limit of a random matrix we find 
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that the kinetic equation is determined essentially 
completely by the correlator G{X ~Y), which de- 
fines correlation between properties of the Hamil- 
tonian describing the "rest of the Universe" at 
two different locations of the simple subsystem. 
Besides the thermodynamic temperature of the 
"rest of the Universe" one needs essentially only 
two more quantities to define the character of the 
dissipative dynamics, the "intensity" of the cou- 
pling controlled by the spreading width (we 
avoid using the term strength of the coupling for 
reasons discussed in Section 2) and by the so 
called correlation length Xq, which defines the 
typical momentum exchanged between the two 
subsystems. In the usual ?i ^ limit, the friction 
and diffusion coefhcients depend explicitly on ?i, 
see Eqs. which seems inconsistent. The 

resolution of this apparent puzzle is simple, one 
cannot simply consider the limit ?i — > 0, while 
keeping Xq and fixed; the corresponding com- 
binations entering in these relations could be held 
fixed, in particular T^II/Xq. 

As we have seen, the freedom of choice our ap- 
proach gives us by allowing various shapes for 
the correlator G{X — Y), leads to new interest- 
ing cases, such of those of Levy flights dynamics 
and turbulent diffusion, and perhaps some other 
cases, which at the moment we do not suspect 
yet. 

Even though the limit of a wide band 
{Kop{e) 3> 1) could prove more then adequate 
for many practical applications, it would be sat- 
isfying to attempt the calculation of the influence 
functional using alternative techniques, like the 
super-symmetry method [|l5|, which do not rely 
on the 1 /TV-expansion. We did mention in previ- 
ous publications that the present approach seems 
amenable to describe localization as well, when 
the finite bandwidth of the random matrix is ex- 
plicitly taken into account. In that case the in- 
fluence functional emerging from this model has 
the structure discussed in Ref. p^ . 

It would also be interesting to determine under 
what general requirements one can derive such 
kinetic equations, in particular by not using ran- 
dom matrix models at all. In spite of the fact 
that random matrix models have proved by now 
to be "universal" , we have seen that to some ex- 



tent even they seem to be too specific. 
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